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Abstract. DM-DE coupling can be a phenomenological indication of a common 
origin of the dark cosmic components. In this work we outline a new constraint to 
coupled-DE models: the coupling can partially or totally suppress the Meszaros effect, 
yielding transfered spectra with quite a soft bending above kj lor eq . Models affected by 
this anomaly do not show major variation in the CMB anisotropy spectrum and it is 
herefore hard to reconcile them with both CMB and deep sample data, through the 
\ same value of the primeval spectral index. 
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1. Introduction 

There can be little doubts that a tenable cosmological model must include at least two 
dark components, cold Dark Matter (DM) and Dark Energy (DE); yet only hypotheses 
on their nature exist, most of them assuming that DM and DE are physically unrelated 
and that their similar densities, in today's world and just in it, are purely accidental. 

Attempts to overcome this conceptual deadlock were made by several authors 
suggesting, first of all, that DE has a dynamical nature pQ (for a review see [2] and 
references therein). An alternative idea is that DE is a phenomenological consequence 
of the emergence of nonlinearity; this appealing option was repeatedly considered (see, 
e.g., [3] and references therein), but is far from being shown and leaves however apart 
the question of DM nature. Interactions between DM and dynamical DE [I] (see also 
[5]) might partially cure the problem, keeping close values for their densities up to large 
redshift. This option could also be read as an approach to a deeper reality, whose 
physical features could emerge from phenomenological limits to coupling strength and 
shape. 

A longer step forward was attempted by [6], suggesting that DM and DE derive from 
a single complex scalar field, being its quantized phase and modulus, respectively. The 
complex field could be the one responsible for CP conservation in strong interactions, 
within a scheme similar to Peccei & Quinn framework [7j (see also [8]). At variance from 
previous suggestions, which introduce parameters and aim at limiting them through data 
fitting, this option - dubbed dual-axion model - cuts the available degrees of freedom, 
including as many parameters as a standard-CDM approach. It is then quite appealing 
that its reduced parameter budget is sufficient to fit quite a number of observational 
constraints [9], still allowing for a common nature of DM and DE and for a specific 
shape of interaction between them. 

In this paper, however, we keep on the phenomenological side and discuss generic 
constraints to DM-DE interactions. This discussion will have a fallout also on the 
dual-axion approach, which does face a problem, because of the feature of the DM-DE 
coupling it causes. 

A coupling of baryons with DE is ruled out by observational consequences similar 
to modifying gravity. Limits are looser for DM-DE coupling, whose consequences can 
be appreciated only over cosmological distances. It should also be outlined that forces 
acting within the dark sector could modify predictions on high concentration DM lumps. 
There can be little doubts that cold DM particles, feeling gravity only, give them NFW 
profiles. Yet, observational data do not lend much support to this shape, for any scale 
range, and direct interaction between DM particles is severely constrained also by recent 
data. 

This is a further reason to consider DM-DE interactions, either as a fundamental 
theory or as an effective framework to approach deeper physics. It is then important 
to devise any observational limit to such interactions and, in this paper, we outline 
further constraints to its shape; they are consequences of the early behavior of density 
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fluctuations, over scales destined to evolve into non-linear structures. 

Fluctuations over such scales enter the horizon before matter-radiation equality 
and their growth is initially inhibited by the overwhelming density of the radiative 
component, then still behaving as a single fluid together with baryons. While 
fluctuations in the fluid behave as sonic waves, self gravitation of DM is just a minor 
dynamical effect is respect to cosmic expansion. This freezing of fluctuation amplitudes 
until equality is known as Meszaros effect. 

The main point we wish to outline here is that DM-DE coupling can damp Meszaros 
effect, so that the rate of fluctuation growth, between the entry in the horizon and 
equality, is significantly enhanced. As a matter of fact, fluctuation freezing is essential, 
in shaping the transfered spectrum, which peaks on the scale kh or ,eq entering the horizon 
at equality. At smaller mass scales (k > kh or ,eq) the spectrum declines because of the 
increasing duration of the freeze. 

The freezing or its damping have modest consequences on the evolution of 
fluctuations in baryons and radiation, evolving then in the sonic regime. What we 
shall therefore find are significant changes in the transfer funcions, while CMB spectra 
keep almost unaffected. 

Constraints to coupled DE models arise from both linear and non-linear effects. 
It has been known since long that coupling may cause a 0-MD epoch after matter- 
radiation equality (see, e.g., pZ3]). This changes the (comoving) distance of the last 
scattering band. In order to fit data, the present value of the Hubble parameter H a 
needs then to be increased. Limits on H Q turn then into limits to the coupling. 

Limits to the coupling, in the case of a Ratra Peebles [11] self-interaction potential, 
where also found in [12], by studying halo concentration distibution. 

The feature outlined in this work affects the transfer function, leaving almost 
unaffected CMB anisotropies. Discussing how the transfer function is affected by DM- 
DE coupling is the main aim of this technical paper. We shall also exhibit CMB angular 
spectra, to confirm that they suffer just marginal changes. 

No general data fitting, constraining parameters and/or showing specific model 
advantages, will be made here. In fact, what we wish to outline is a major effect, 
which allows to discard a class of models, a priori. This is why we keep to cosmological 
parameter values ensuing from WMAP3 best-fit [IB], although deduced by assuming 
a ACDM model. In particular, we shall take an overall density parameter Q = 1; 
the present value of the cold DM (baryon) density parameter will be f2 0jC = 0.224 
(fl ,b — 0.044); the dimensionless Hubble parameter will be h — 0.704; the primeval 
spectral index, when not taken as a free parameter, will be n = 0.947. 

Within this frame we shall consider a self-interacting scalar field, causing cosmic 
acceleration when its pressure/density ratio w = PdeI Pde falls in the range (—1,-1/3). 
Quite in general, it is 

Pde = Pk,DE + P Pi de = <fi 2 /2a 2 + V ((/)), Pde = Pk,DE ~ P p ,de , (1) 
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so that it is— 1/3 ^ w > — 1 when dynamical equations yield Pk,DE/V "C 1/2. Here 

<is 2 = g iiv dx ii dx v = a 2 (r)(-dT 2 + dxidx 1 ) (i = 1, .., 3) (2) 

is the background metric and dots indicate differentiation with respect to r (conformal 
time). The w ratio exhibits a time dependence set by the shape of V(<ft). Much work 
has been done on dynamical DE (see, e.g., [2] and references therein), also aiming at 
restricting the range of acceptable u>(r)'s, so gaining an observational insight onto the 
physics responsible for the potential V(4>). 

Our analysis here will however be restricted to SUGRA potentials [14] 

V((f>) = (A a+A /r) exp(47r0 2 /^) (3) 
admitting tracker solutions. Here m p = G -1 / 2 is the Planck mass. This will enable us 
to focus on peculiarities caused by the coupling. Let us also remind that, once the DE 
density parameters Qde is assigned, either a or the energy scale A, in the potentials 
([3]), can still be freely chosen. In this paper we show results for A = 10 2 GeV. 

The SUGRA potential, at least in the absence of coupling, yields an excellent fit of 
observational data [T5] . We tested the effects of coupling for a number of values of the 
scale A, from 10 to 10 4 , and also changing the shape of the potential into Ratra-Peebles. 

In the former case we find just marginal shifts. In the latter one, quantitative 
changes can be significant. The overall behavior is however identical and this potential 
is known to yield a poor fit to CMB data, unless quite a small A scale is taken, so 
spoiling its physical appeal. 

Within this frame we aim at focusing problems and showing the quantitative 
consequences of different options. 

2. Dynamical equations 

Let us then start from the background equations for DE and DM when the metric is 
(T2]), reading 

+ 2-0 + a 2 V^ = C{(f))a 2 p c , p c + 3-p c = -C(0)0p c (4) 
a a 

where we set 




A possible time dependence of the coupling strength was considered but not deepened 
since the early work of [T5]. The dual-axion model naturally predicts a coupling 
C = 1/0, consistent with eq. ([5]), if e = — 1 and (3 ~ 0.244. Quite in general, for 
dimensional reasons, C can be expressed through products of and 6 with a+b = — 1. 
During most cosmic evolution is a monotonically increasing function, so that the sign 
of the exponent of tells us whether the coupling was stronger of weaker in the past. 
(An exception can be recent times, as the exponential term in the SUGRA potential 
can cause a re-bounce of when it approaches m p ; our arguments here concern much 
earlier times, when it is safely < m p ). 
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An expression of C, made by a polynomial including terms with different powers 
of </>, could select a peculiar epoch to have then a weaker or stronger coupling. Such an 
option, however, is clearly ad-hoc and does not seem to deserve further investigation. 

An explicit dependence of C upon time would be hard to reconcile with the Lorentz 
invariance of the Lagrangian mass term 

C c = -B(cj))m x xx (6) 

setting the coupling between the DE scalar field <p and a spinor field x supposed to 
yield DM, as C(4>) = d(\nB)/d(j). (Notice that B((p)m x is the time-dependent mass 
of DM quanta). A dependence of C on seems therefore the only way to instaure a 
time-dependent coupling. 

Let us then describe fluctuation equations. In the period when Meszaros effect 
occurs, DM, photons and baryons can be treated as fluids; (massless) neutrinos, instead, 
are not a fluid. Fluctuations in a generic fluid with p/p = w and 5p/5p = c 2 s fulfill the 
equations 

5 = -(1 + w)(kv + h/2) - 3(c 2 s -w)(H- C<p)5 - (1 - 3w)(C<p + C'jxp) 

v = -(l- 3w)(H - C6)v + — —kS v -ha- kC cp . (7) 

v /v ' 1 + w l + w 1 + w 

Here 5 = Sp/p, v = iuik l /k, (p +p)c = —(hi ■ kj — 5y)(T u — 5ijT^/3) (u{ are the space 
components of the velocity field in the fluid and is its stress-energy tensor) and 
TC = a j a, while the DE field 

0(r,x) = o (r) + <p(r,x) (8) 

is split into a background component (f) , coinciding with the <fi field obeying the 
eq. d3|, and the space-dependent fluctuation if] in eqs. ([7]), Fourier components of ip are 
considered which fulfill the equation: 

(p + 2H<p + k 2 <p + a 2 V£<p + <j)h/2 = Ca 2 p c 5 c + Cy Pc <p , (9) 

while fluctuation self-gravity is fully accounted by h, obtained by integrating the 
equation 

h + Hh= -8nG {a 2 p[n ib (l + 3c 2 s )5 lb + tt c 5 c ] - 2a 2 V^(p + (pip) (10) 

In the absence of coupling (C = 0), we can take v = in eqs. (|7j) without loss of 
precision, and face the dynamical problem though a single first order equation. This is 
no longer true when DM particles can be pushed by DE forces. Then the DM equation, 
on scales below the horizon and on times before matter-radiation equality, reads 



a 4 7T ~, 



a m p 



5 c -47rGa 2 p{ 



1 + ^ 2 (0) 



l + 3c 2 



, (11) 



Here, p is the overall density; fl c and f2 7 h are time dependent density parameters for 
DM and the baryon-photon fluid. To show eq. (fTTj) one needs a little algebra, which 
will be postponed to the end of the section. 
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Still from eqs. fl7|), we can also work out the equations for the photon-baryon fluid, 
by taking C = and 

c 2 = 5 Pj /(5 Pb + 5 Pj ) = [3(1 + 306/4O,)]- 1 , (12) 
w 7 fe(a) = P 7 /(Pb + P 7 ) = [3(1 + Qb/Q^)]' 1 . (13) 

Here and f2 7 are the time dependent baryon and photon density parameters. The 
baryon component can be responsible for a shift of from 1/3; although initially small, 
it can approach 20-25 % at the eve of recombination. Notice then that a non-vanishing 
factor 1 — 3w ^ Q b /Q ry keeps a direct influence of ip on baryon-photon fluctuations. 

This set of equations enables the reader to build a simplified numerical algorithm, 
directly testing the suppression of Meszaros' effect. 

Clearly, to study the later evolution, since the eve of baryon-photon decoupling 
a full kinetic treatment of the radiation is needed. This will be used to confirm that 
CMB anisotropics are just marginally affected, but is unessential to focus on the greater 
changes occurring to the transfer function. 

Let us now summarize the procedure to obtain eq. (jTTj) ; a reader unintersted in it 
can skip the rest of this section. 

The starting point are again the eqs. ([7j). Together with them, let us consider again 
eq. (EJ). There, any mass-like term multiplying 0, in comparison with k 2 , is negligible; 
then, before equality, it yields 

(p + (2/r)<p + k 2 (p = Ca 2 p c 5 c - <ph/2 (14) 

Here, two kinds of time dependence must be compared, over fluctuation and Hubble 
time scales. Accordingly, we can express (p as sum of rapidly and slowly varying terms 
by actually summing up the (rapidly varying) general integral of the equation obtainable 
by equating to zero the l.h.s. and a (slowly varying) integral obtained by equating the 
last term at the l.h.s. with the r.h.s. . If we then time-average over fluctuation time 
scales, only the latter contribution survives and 

{(p) ~ p (Ca 2 p c 6 c - 4>h/2) (15) 

Let us recall that dropping the contribution of fluctuating terms is the standard 
procedure to obtain an analytical description of Meszaros' effect. The point here is 
that, while all ip derivatives can be dropped, there is a slowly varying contribution to ip 
which cannot be soon disregarded. 

However, if the expression <HM is used to replace ip in eqs. (|7|) for CDM, yielding 



5 c + kv c + ^ = -2(3 2 ^- 
2 dr 



Jkrf 



+A ^l 2 i{6^) (16) 



it becomes clear that also the slowly varying ip contributions can be dropped. In fact, 
the first term at the r.h.s. contains a division by (kr) 2 , which is the squared ratio 
between long and short timescales and, altogether, it is ~ 0[5/(kr) 2 ] <^ 5. We must 
then acknowledge that a time derivative of a quantity Q, varying over the Hubble time 
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scale, is ~ 0{Q/r). The second term at the r.h.s. is then ~ (9[(0/m p ) 1+<E /i/(/cr) 2 ], while 
we took 1 + e > and, in the epoch considered, <C 1. It must then be even smaller 

than the the first one. 

Setting then to zero the l.h.s. of eq. (flEI) . we obtain the relations 



5 C + kv c + 



h , h 

- ~ , kv c ~ -5 C - - . 

In turn, the second eq. (J7J), using this latter equality, yields the relation 

kv v = (Cj> - H) \-8 c - h/2] - k 2 C^ , 



(17) 



which can be replaced in the former eq. (flTj) . together with eq. (TlOT) and ( fT5l) . obtaining 

H - C<j>] 5 C - ArtG [a 2 p[n, /b (l + 3c 2 s )5 lb + Q C 5 C ] - 2a 2 Vjy> + - C 2 a 2 p c 5 c = 

Neglecting here if fluctuations as source terms - DE yields a minor contribution to the 
overall density, as shown in Fig. [6]- and using eq. (jSJ), eq. (1TTI) is soon obtainable. 



3. Overcoming the freeze 

Let us then focus our attention on eq. fill I) and outline first that, in average, the term 
^76^76 (1 + 3c 2 ) almost vanishes, as sonic waves fluctuate. Then, in the absence of 
coupling (/3 = 0), the self-gravitation term oc 5 C is also damped by ft c <C 1- If the 
photon-baryon term is then neglected, the increasing mode (approximately) reads 

5 C oc 1 + 3y/2 (19) 

with y = 3w 1 bO J /a eq (see, e.g., [IE]) and this yields a growth from horizon to equality 
never exceeding a factor 2.5 , that we shall approximate as 5 C oc a 1 / 4 , according to 
numerical outputs. Meanwhile, above the horizon, in a synchronous gauge, 5 C oc a 2 . 
Hence, for k > kh or ,eq, the growth is slowed down by a factor oc a^ 4 (A;), while the scale 
factor when k passes through the horizon, ah(k) oc k . Altogether, for k > kh or ,eq, 
we expect a transfered spectrum P(k) = Ak n T 2 (k) ~ Ak n ~ 3 ' 5 (T(k) is the transfer 
function). 

It is then easy to see what can change because of the coupling. The coefficient 
of the friction term oc 5 C is certainly reduced and can even invert its sign. Even more 
significantly, the term [l + 4/3 2 (0) /3]Q C can attain or overcome unity, not only because of 
the greater size of $(<p), but also thanks to the modified time dependence of Q c . It turns 
out, in fact, that a decreasing dependence of $ on (ft (e < 0) yields higher Q c values in 
the relevant redshift interval. In general, is smaller at earlier times and e < causes 
a stronger coupling in the past. 

Let us however debate first the constant coupling case (e = 0). Available data 
set then a constraint (3 < 0.1-0.2 [10] . [12] (see also [T7]; beware of the different 
coupling definition), limiting the acceptable discrepancy of DM and DE evolution, after 
recombination, from uncoupled models. The most direct effect, in this case, concerns 
large scales entering the horizon late, when such discrepancies occur. The low-/ plateau 
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Figure 1. Scale dependence of the density parameters of the various components in 
coupled DE models with constant coupling. This plot shows also the displacement of 
z eq and, henceforth, of kh or ,eq as j3 increases: thicker (thinner) curves refer to (3 = 0.01 
(0.0244). 




Figure 2. Best fits of SDSS data for constant (3 from (solid line) to 0.25 (dotted 
line). Different lines correspond to a (3 increase by 0.05 . The vertical dotted line yields 
the scale of Cio- 
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Figure 4. n intervals for increasing (constant) coupling strength 

of the CMB anisotropy spectrum can then undergo a modified ISW effect, while some 
changes in the C\ behavior, up to the first peak at I ~ 200, can be compensated by 
slightly modifying n, h and other parameters. 

But the most significant shifts occur on the transfer function. Its slope, at k > kh or ,eq 
and up to a scale k ~ 0.1/ZiMpc" 1 , where non-linearity effects become important, is 
slightly distorted as (3 increases. The main effect, however, is a progressive displacement 
of kh or eq itself. In Figured] we compare the Q evolution in models with different constant 
coupling, showing the significant displacement of the crossing between Q r and Q c , when 
different /?'s are taken. 

We can show the impact of this displacement by fitting transfered spectra, over 
these scales, with the Luminous Red Galaxies sample data from the Sloan Digital Sky 
Survey (SDSS) [18] and allowing the primeval spectral index n, assumed to be constant, 
to act free parameter. 
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In Figure [2] we show the result of this fit. Transferee! spectra, when kh or ,eq vary, easily 
accommodate deep sample data in the linear range, as the kh orm shift is compensated 
by a slightly smoother slope. All that however requires a non negligible decrease of n, 
and transfered spectra risk to become too high in the region where they should fit CMB 
data (the scale of the 10-pole is indicated in the Figure). 

An attempt to balance this spectral distortion can be made by varying other 
parameters. It is then significant to consider Figure [31 showing deep sample data 
vs. transfered spectra computed with n's within 1- and 2-cr's from the best-fit. In 
the Figure we took (3 = \J?>/ix /A ~ 0.244. Then, in Figure HI we show the 1- and 2-cr 
range of n, for constant coupling strength. Let us also recall that likelihood analysis 
showed that models with constant (3 <~ 0.2 do not exhibit severe disagreements with 
data. 

The effect on CMB data fitting arising from n values distant from unity can be also 
directly inspected in Figure [5j According to it, we can examine Fig. @] assuming to be 
viable only those models which, at the 1—a level, admit n >~ 0.85 . With reference to 
this admittedly qualitative criterion, we shall now consider the variable coupling case. 

First of all, when coupling varies, the evolution of the density parameters exhibits 
significant discrepancies from the constant coupling behavior. The point is that, at 
high redshift, they further strengthen DM self-gravity, coherently with higher (3 effects. 
Results of a numerical integrations, illustrating this issue, are shown in Figure El 
where we compare the redshift dependence of density parameters in constant (e = 0) 
and strongly variable (e = —1) coupling models, allowing to appreciate a substantial 
enhancement of Q c at the eve of equality. Altogether [1 + 4/? 2 (0)/3]f2 c is significantly 
greater, and the freeze of 5 C , between horizon entry and equality, is almost canceled. 

This explains the behaviors shown in Figure 0, which exhibit one of the main 
findings of this work. These plots are obtained from a numerical evaluations of 5 C , in 
models with different (3 for e = or -1. In the former case, the high-2; effect of coupling 
is marginal. In the latter one, the fluctuation growth is substantially enhanced, more 
significantly for greater k values. 

Accordingly, we can expect modifications in the transfer function and in the fit 
of observational data. Our general conclusion is that a time dependence of DM-DE 
coupling, making it stronger at higher redshift, can prevent DM fluctuations to have a 
stationarity period after their entry in the particle horizon, so causing large modification 
of the transfer function. 

4. Results 

These expectations will be tested by using an algorithm solving the whole set of 
dynamical equations. To this aim we can use our extension of the public programs 
CAMB, or our own code, with identical outputs. 

Results for transfer functions are shown in Fig. [8] for (3 = 0.1 and 0.244, and a 
variety of values of e. 
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Figure 5. Anisotropy spectrum of the ACDM model yielding the best fit to WMAP3 
data compared with the spectra for coupled models with /? = 0.1 and e = 0, for n = 1 
and n = 0.7. Already in the former case some difference exists, but no major qualitative 
changes occur; by adjusting other model parameters one can expect a reasonable fit to 
data. Taking n — 0.7, any fitting to CMB anisotropy data is apparently excluded. 



The suppression of fluctuation freezing is obviously stronger for greater (3 (and 
increasingly negative e values). For e = — 1, enclosing the case C = 1/0 when {3 = 0.244, 
the steepness of the transfer function, for k > kh or ,eg is much reduced. The effect is still 
significant also for e = —0.5, namely when j3 = 0.244. 

A further effect shown by these plots is a significant displacement of the scale 
where T begins its gradual descent. As a consequence, different coupling laws may 
cause displacements on the scale where transfered spectra peak. 

Notice that, while this occur, the CMB anisotropy spectum keeps quite a reasonable 
behavior, if n ~ 1, as is shown in FigureEl similar to Fig.[5]but for variable coupling; here 
we compare the WMAP3 data on the anisotropy spectrum with the spectra obtained 
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Figure 6. Rcdshift dependence of density parameters fl in coupled DE models with 
constant and variable coupling. The two panels refer to different values of /3. In both 
panels we show fi's for e = (C oc f3/m p , thinner lines) and e = — 1 (C oc [3/cf), thicker 
lines). Figures are rather intricate and their reading may begin from dashed lines, 
yielding Sl r = il 7 + £1„, which arc almost independent from the coupling law. A more 
relevant effect occurs on f2 c (dotted lines), whose values, in the case of variable coupling, 
exceed those of constant coupling until equality. The most relevant effect occurs for 
DE (solid lines), whose contribution to the overall density is enhanced by several order 
of magnitude by variable coupling. DE fluctuations will fade after the entry in the 
horizon; however, high SIde values increase their impact on other components before 
their disappearance. 
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Figure 7. Evolution of DM fluctuations in a time interval enclosing the entry in the 
horizon and the matter-radiation equality. In all cases, in the presence of DM-DE 
coupling some modification occurs. They are however almost negligible for constant 
coupling, while, for coupling oc <^> _1 , Meszaros' freezing is almost completely canceled. 
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Figure 8. Transfer functions for different behaviors of DM-DE coupling with redshift 
and/or for different coupling normalization. The case e = corresponds to redshift 
independent coupling intensity. The case e = — 1 with (3 = 0.244 correspond to a 
coupling C = \/4> . Besides of the different slopes, notice the dependence on the model 
of the bending scale and, in particular, its dependence on the coupling strength, also 
in constant coupling models (dash-dotted lines). 

for (5 = 0.1, e = — 1 and two n values: n = 1 apparently allowing a reasonable fit, 
and n = 0.7, as needed to obtain a reasonable transfered spectrum. Figure in fact, 
is a direct evidence that, in order to recover a fair slope of the transfered spectrum in 
the scale range where structures accumulate, small primeval n values are unavoidably 
required. In order to perform an evaluation of the effect, we actually built transfered 
spectra and compared them again with SDSS data. In this way we find n ~ 0.5-0.7. 
Fig. [9] shows that this spoils the fit with G\ data. 

Spectra with ordinary downward bending fit both deep sample and CMB data with 
n ~ 1. On the contrary, spectra with low n's and standard oVs cause greater C\ still in 
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Figure 9. Anisotropy spectrum of the ACDM model yielding the best fit to WMAP3 
data compared with the spectra for coupled models with f3 = 0.1 and e = 1, for n = 1 
and n = 0.7. As in Fig. 03 in the former case one can expect to recover a reasonable fit 
to data by adjusting other model parameters . Taking n = 0.7, a value just acceptable 
to fit deep sample data, any fitting to CMB anisotropy data is apparently excluded. 



the Sachs & Wolfe plateau, while reducing the relative height of the C\ peaks (see again 
Fig. [91 These effects could be partially compensated by an adjustment of other model 
parameters, whose search is out of the scopes of this work. Finding n values below 
0.8-0.9, however, clearly means that we are dealing with unlikely physical frameworks, 
hardly allowing to fit CMB and deep sample data simultaneously. 

In the case e ^ 0, the discrepancy from unity of the spectral index n, assumed to be 
constant, is mostly a measure of the distortion caused by the suppression of Meszaros 
effect, which overwhelms the effects of the displacement of kh or ,eq, already considered 
in the (3 = const, case. Moreover, such discrepancy is a significant estimate of the 
distance of the model from uncoupled physics. Our fits, shown in Figures [10] and [Til 
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Figure 10. Model comparison with SDSS digital survey data. Different curves refer 
to different values of e (as in previous Figure) with the solid line (e = f ) essentially 
coinciding with an uncoupled model. Constant coupling models (e = 0) are described 
by dot-dashed curves. Negative e's yield a further decrease of n. The vertical dotted 
line is the approximate scale where the Sachs & Wolfe C\ plateau begins. Constant 
coupling causes a rise of Cio by a factor ~ 1.8. A further factor ~ 2 arises from a 
coupling C — 




0.001 0.01 0.1 
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Figure 11. As fig. [TUl for a smaller coupling intensity. The effect on the Cio scale 
is much reduced, for constant coupling; however, variable coupling yields a dramatic 
decrease of n and, for e = —1, the level of (3 — 0.244 is almost attained. 

are complemented by 1- and 2-cr intervals around best-fit n values (Figure [T2"j) ; they 
are an indication of which models, by adjusting other parameters, might be susceptible 
to approach the observational scenario. The behavior of transfered spectra in respect 
to data, when taking n values at 1- and 2-cr from best fits is also shown in Figure [13] 
(similar to Fig. [3]) . 

If we take again n = 0.85 at I— a as a threshold to discard a model, no e < coupled 
model is allowed with f3 = 0.244, while e < —0.16 are also inhibited with /3 — 0.1 . At 
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Figure 12. 1- and 2-cr intervals of n, when e varies, for /3 values. 
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Figure 13. If values of n at 1- or 2-cr's from best fits are taken, spectra are significantly 
modified. Here we show the effect in the case with C = 0.4( 7 r/3) 1 / 2 /y^. 
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2-cr's the situation is not much improved for (3 = 0.244, while lower values of e are 
admitted for (3 — 0.1 . In particular, a model with C = 1/0, as the dual-axion model, 
lays outside of the range indicated. 

The analysis was extended here to models with positive e, for which coupling rises 
while (f) increases. A large deal of these models is apparently allowed. 

5. Conclusions 

The quest for models fitting observational data and avoiding fine tuning and coincidence 
suggests to test cosmologies where DM is coupled to dynamical DE. It is then important 
to recognize that quite a few models, with constant or variable DM-DE coupling, are 
consistent with observational data, although their likelihood might be slightly smaller 
than uncoupled models. 

Constraints on these models were known to arise from a number of linear and non- 
linear effects. The main linear anomaly is the existence of a prolongated period, after 
matter radiation equality, when the </>-field energy affects the expansion rate (the <p- 
MD epoch). As a consequence, the comoving distance of the last scattering band can 
be different from ordinary models. In principle, this can be compensated by varying 
the present value of the Hubble parameter H . There are however severe limits on H , 
which turn into limits on the coupling strength, discussed in [9] . The same feature causes 
also a displacement of the wave number kh or ,eq, corresponding to the scale which enters 
the horizon at matter-radiation equality. Here we discussed also this effect, which was 
however expected. 

Constraints on coupled DE models were also found by studying their non-linear 
evolution. In [32] it is shown that coupling can affect the halo concentration. Although 
this effect could be softened by using suitable self-interaction potentials, it yields a limit 
(3 <~ 0.15-0.20, for the coupling, when the self-interaction potential is Ratra-Peebles 

M. 

In this paper we outlined a new effect, that DM-DE coupling may cause in a 
class of models. A prolongated period, between the horizon entry and the equivalence, 
when DM fluctuation growth stagnates, is essential in shaping transfered spectra P tr (k). 
Because of this Meszaros's effect, Pt r {k) bends at k > kh orfiq . In this paper we outline 
that this stagnation period can be partially or totally suppressed by the coupling of 
DM with the <fi field. As a result, transfered spectra exhibit a much softer bending at 
k > khor,eq- We also outline that, while this occurs, the dynamics of sonic waves in the 
baryon-photon fluid is only marginally affected. Accordingly, while transfered spectra 
suffer major changes, CMB anisotropies are almost invariant. 

Using a single primeval spectral index n, we can fit both CMB and deep sample 
data, in the presence of standard Meszaros effect. In the presence of the above softening, 
the values of n required to fit CMB and deep sample data become badly discrepant. 

In this work we discuss which class of coupled DE models are affected by this 
anomaly. If the coupling is constant or proportional to a positive power of 0, the 
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anomaly is absent. On the contrary, models where the coupling is proportional to an 
inverse power of <ft are at risk. 

We gauge the impact of this anomaly by evaluating the value of n required to fit 
just deep sample data with each model. When n <~ 0.85, it is legitimate to believe that 
the model likelihood is highly suppressed. The search of such n value should therefore 
be preliminary to any attempt to reconcile a coupled model to the whole set of CMB, 
deep sample (and other) data. 

Admittedly, apart of conceptual reasons, within the present observational 
framework there lacks any specific phenomenological push to invoking a DM-DE 
coupling. Yet, uncoupled models, even apart of their conceptual weakness, cause a 
number of questionable predictions, e.g. NFW profiles. Interactions within the dark 
side were often advocated to cure such difficulties. An impact on profiles was actually 
shown to exist, but new form of coupling need to be inspected. 

Furthermore, fresh data on the redshift dependence of p c , pb and Pde, at z ~ 1-5, 
might soon be available, if experiments like DUNE |20j will become operational. The 
discovery of an anomalous scaling of p c , for instance, would set a strong prior, completely 
biasing likelihood distributions, just as priors on the value of h (Hubble parameter) 
suppress the likelihood of SCDM models with h ~ 0.4 , which would otherwise allow a 
reasonable fit of large sets of data [2TJ . In particular, models with rising coupling (e > 0) 
could become an important option. Such coupling behavior could directly arise from 
suitable microphysics or be a phenomenological description of a complex underlying 
physics. 
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